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Resume 

The aim of this paper is three-fold. First to estabhsh in a clear and rigorous 
way a formula proposed heuristically by Mehlig- Wilkinson for the Metaplec- 
tic Operators corresponding to a given Symplectic Transformation in classical 
Phase-Space. 

Second this formula is applied to the study of Quantum Recurrences, which 
has attacted great interest in recent years (see [SEl for a complete account of 
the recent approaches). The return probability is given by the squared mod- 
ulus of the overlap between a given initial wavepacket and the corresponding 
evolved one ; quantum recurrences in time can be observed if this overlap is 
unity. We provide some conditions under which this is semiclassically achieved 
taking as initial wavepacket a coherent state localized on a closed orbit of the 
corresponding classical motion. 

Third, we start a rigorous approach of the "quantum fidelity" (or Loschmidt 
Echo) : it is the squared modulus of the overlap of an evolved quantum 
state with the same state evolved by a slightly perturbed Hamiltonian. It 
has attracted a great interest in the last decade, for the purpose of "Quantum 
Chaos" problems, and to quantum computation analysis. However the results 
are most of the time not entirely conclusive, and even sometimes contradic- 
tory. Thus it is useful to start a rigorous approach of this object. 
The decrease in time of the Quantum Fidelity measures the sensitivity of 
quantum evolution with respect to small perturbations. Starting with suitable 
initial quantum states, we develop a semiclassical estimate of this quantum 
fidelity in the Linear Response framework (appropriate for the small pertur- 
bation regime), assuming some ergodicity conditions on the corresponding 
classical motion. 
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1 Introduction 

The revival of quantum wavepackets along the time evolution has attracted much 
recent interest (see [36] and references therein contained) , and we want to stress here 
how much can be demonstrated about this topic in the semiclassical limit, starting 
from a wavepacket (coherent state) conveniently localized around a point on a closed 
classical orbit (Section 3) It consists in considering the overlap of a time evolved 
quantum wavepacket with its initial state at conveniently chosen times called "re- 
vival times". 

In addition a large physical literature has reported (essentially heuristic) re- 
sults about the so-called "quantum fidelity" (also called Loschmidt echo) in the 
perturbative and/or semiclassical limit (see P0,[21, U\, |H| ; [IH] ; [IH] ; [H] ; [T^ - 

EDI, EH, 122, [23,1211, 12ni,|SIl,E21,E3,ElI, ESl,EIl,EHlEni,ffl,|iIl, 1121,113, 

111,113, |ini,|lI|,|lHl,|13,|Sni,ISIl)- The idea, which goes back to Peres [SH is to 
study the behavior in time of the sensitivity of quantum evolutions with respect to 
small perturbations of the Hamiltonian, mainly in situations where the underlying 
classical motion is "chaotic" . 

More precisely, one considers the overlap of the perturbed quantum evolution of 
some given initial wavepackets with that under the unperturbed quantum dynam- 
ics, in particular the decay properties in time of this overlap. 

Here we shall consider exact results about this behavior in the "linear response" 
framework, starting with a bunch of initial wavepackets taken as eigenstates of the 
unperturbed Hamiltonian, assuming some ergodicity or mixing assumptions about 
the underlying classical flow, and taking the semiclassical limit h —>■ 0, (Section 3). 

We note that these results constitute a first step towards a completely rigorous 
understanding of these topics in the semiclassical regime. 

The paper is organized as follows. In section 2, we give a proof of a beautiful 
formula proposed by Mehlig- Wilkinson [23 for the metaplectic operators associated 
to a given symplectic map. (This topic is being developed and completed in a recent 
work in collaboration with D. Robert). This formula appears to be very useful in the 
semiclassical study of the Quantum Revivals that we develop in Section 3. Section 

4 contains some results about the revivals of coherent states in the case of time- 
periodic Hamiltonians that are quadratic in coordinates and momenta. In Section 

5 we give some result about the Quantum Fidelity in the semiclassical regime, in the 
Linear- Response framework, and recover a result derived heuristically by Prosen and 
Znidaric |23- Section 6 is an appendix containing intermediate lemmas necessary 
for the proof of Theorem 3.6. In Section 7 we give some concluding remarks. 
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2 A proof of the Mehlig- Wilkinson Formula 

In a recent paper, Mehlig and Wilkinson proposed (heuristically) a beautiful 
formula expressing the metaplectic operators corresponding to suitable symplectic 
matrices M not having 1 as eigenvalue. In a paper in preparation with D. 

Robert (Pl]), we give a developed and generalized proof of this formula, including 
the calculus of the phase, and the case where the Symplectic Matrix has 1 as eigen- 
value. Here we present a simple proof, closer to the physical intuition, in the same 
spirit as in |14 but without expliciting the phase factor. 

A metaplectic transformation implements in the quantum world the symplec- 
tic transformations Sp{2n) (canonical transformations) of classical mechanics. To 
a given a symplectic transformation M, it corresponds unitary operators R{M) in 
the Hilbert space of quantum states acting covariantly on the canonical operators 
of position and momentum Q, P in the following way : 

R{Mr 9 ^ R{M) = M ( 9 ^ (2.1) 

plus some algebraic conditions ensuring that M i— R{M) defines a projective repre- 
sentation of Sp{2n), with sign indetermination. For any M G Sp{2n) we can find a 
C^-smooth curve Ft, t G [0, 1] in Sp{2n) such that Fq = 1 and Fi = M. Ft is clearly 
the linear flow defined by the quadratic Hamiltonian Ht{z) := ^z.Stz with 

St = -JFtFt (2.2) 

where 

If Ut is the quantum propagator associated with Ht, we define 

R{M) = Ui 

which obeys equ. (j2.H ). If two sympectic paths Ft and join 1 to M, then we 
have Ui = ±U[. 

Definition 2.1 Let z := {q,p) be a generic point in classical phase space Z := M^". 
The covariant symbol of an operator A is a function A'^{z) defined in the whole 
phase-space M^" such that 

A = j dzA\z)f{z) (2.4) 
with T{z) being the Weyl Heisenberg operator 

T(.):=exp|'^?^«^| (2.5) 
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Theorem 2.2 Assume that M doesn't have 1 as eigenvalue. Then the covariant 
symbol of R{M) is given by 

^^'^ = |det(ll-M)|V2 ^-P (^^•^-) (2-6) 

where 'Jm is a complex number of modulus 1, and A is defined by 

A:=^{M -t)-\l + M) (2.7) 

and J is the symplectic matrix defined in equ. ( t^./i]) . 
Proof : By the inversion formula we get 

n{X) = Tr(f/if (-X)) (2.8) 

which using the coherent states basis 

:= f{z)^o (2.9) 

where 

:= (7rn)-"/^exp(-a;V2n) (2.10) 

can be rewritten as : 



Q{X) = {2nhr'' J dz{^,, f{-X)U,^,) (2.11) 



In all that follows we shall omit for simplicity the factors 27r and h, that we can 
easily reintroduce in the result. We have : 

f{X)f{z) = e-^-(^'^)/2f (X + z) (2.12) 

with the notation 

a{z, z') = z.Jz' 

Clearly 

Uiifi, = fiz')RiM)ipo (2.13) 

with z' := Mz, so that 

n{X)= /rf^e'"(^'^)/2-'"("''^+^)/2(^y,^(M)^o) (2.14) 



where by Y we denote the variable 

Y ■.= X + z~z' = X + {1-M)z (2.15) 
The scalar product {ify, R{M)ipo) is equal to 

JY.QY/2 
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with Q a matrix that we do not need to make exphcit. Performing the change of 
variable z ^ Y, we rewrite a{z + z',X + z) in terms of X and Y alone. The cross- 
terms vanish, so that we are simply left with 

-a{z + z',X + z) = X.AX + Y.Q'Y (2.16) 

with Q' a matrix that we do not need to make explicit. For the calculus of A, we 
easily check that the quadratic terms in X are : 

-{1+M){l-M)~^X.JM{1-M)-^X = -{l~M)-^X.J{l+M){l-M)-^X (2.17) 

where we have used that, by the symplecticity of M, 

{1 + M)JM = JA4 + J = J{1 + M) (2.18) 

where M denotes the transpose of M. Therefore ()2.15|) is nothing but 

-(1 - My'X.JM{l - M)-'X = -(1 - M)-'X. J^^-^^J _ M)-^X (2.19) 

However, using again the symplecticity of M, we have : 

(1 - M)-\JM - MJ) = J(l + M) 
Therefore (|2.19p is nothing but 

--X.J(11 + M)(ll - M)~^X = X.AX 

Therefore inserting this result into ()2.14j) simply leads to : 

n{X) = fi(0)e*^-^^/2 (2.20) 

where Q{0) is the integral over Y that we do not perform exactly, although it is pure 
Gaussian integral, but we shall content ourselves to determine its absolute value. 

Lemma 2.3 We have 

|n(0)| = |det(]l-M)|-i/2 (2.21) 

Proof : 

(<^„t/i<^,)= / dXU,{X)W^^{X) 



where Ui{X) is the Weyl symbol of the quantum evolution operator Ui, and W^^{X) ■ 
W^,{X-z) = MX-z) with 

M^) ■■= e-'' (2.22) 



being the Wigner function of (fo. Thus 

{ip,,Uiip,) = Ui*%iz) (2.23) 
Let us now calculate / dz\{{pz, f/iv^z)^ by using Plancherel Theorem : 

j dz\{ip,,U,<f,)\' = j dz\{Ur*^o){z)\^ (2.24) 
d;2|^(f/i*$o)(z)P 

where JF denotes the symplectic Fourier Transform. But VL{z) is the symplectic 
Fourier Transform of Ui{z). Thus equ. ()2.24|) becomes 

= j dz\n{z)\''\r<^Q{z)\^ = |l](0)p j dze-'''^ = |r](0)|2 (2.25) 

because A being real symmetric (which can be deduced easily from the symplecticity 

oi M),\n{z)\ = \Q{^)\, VzgM^". 

Let us now calculate the same quantity using the Wigner Transforms : 

f/i^,)p = 2 I du W^t(.)^o(«) ^f(.')ij(Af)^o(«) (2-26) 

= 2 j duW^,{u - z) VrR(M)^.> ~''^^'^J '^''''''P ^'^"^ ~ " ~ ''^'^ 

Now, in order to integrate equ. ()2.26p over z, we perform the change of variable 

z X := z — u 
so that we are left with a simple Gaussian of the form 

2 j dXdu e-^'-(^-^)^ = j du e-^'/2 = | det(l - M-^)\-' (2.27) 

since Y := {M~^ — l)u. But since M is symplectic, it has determinant one and thus 
equ. p.27|) equals 

|det(M- 1)|-^ 

Hence, 

|O(0)p = |det(ll-M)|~^ (2.28) 

□ 

Theorem 2.2 is nothing but the beautiful Mehlig- Wilkinson formula ([SHI) for the 
metaplectic representation. In J3j , we establish it in full generality, including the 
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calculus of the phase factor 7(-F). Assume that F is a symplectic map not having 1 
as eigenvalue. Then : 

where 7(-F) is a complex number of modulus 1 that we do not specify, and A is the 
2n X 2n real symmetric matrix 

A := ^J{F + 1){F - 1)-^ (2.30) 

and J is the symplectic matrix ()2.3|) 
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3 Quantum Revivals 

We consider a quantum wavepachet ip E Ti. = L^(M'^), together with a quantum 
Hamiltonian H selfadjoint in Ti. According to Schrodinger equation, the quantum 
evolution of wavepackets is provided by the unitary group of evolution Unit) '■= 
ex]i{—itH/fl), h being the Planck constant (divided by 27r). The quantum revivals 
of that wavepacket are observed by considering the overlap between the evolved 
wavepacket UH{t)(f with the initial one (p : 

R{t) := {^,UH{t)v>) (3.1) 

Assume now that the initial wavepacket is taken as a Gaussian one, namely a 
coherent state ip^ defined in equ. (2.9-10) for any z := {q,p) E Z := M?"'. 

The coherent states (sometimes written in the Dirac notation \z) instead of (p^) 
have been shown to display remarkable semiclassical propagation properties (see 
|12j ) that we recall here : 

Along a tradition which goes back to Hepp (1974), one can start by "following" 
a coherent state along its semiclassical evolution. We shall establish the following 
result : 

starting from a coherent state \z) at time t=0, its quantum evolution stays, up to a 
phase, close to a "squeezed state" 

centered around the point Zt := (p^iz), ($^ being the classical flow generated by the 
classical symbol H of < H), with a "dispersion" governed by a symplectic matrix Ft 
that we shall make precise later. This approximation is of order 0{h'^) as long as time 
doesn't go beyond the so-called Ehrenfest time Te '■= logh^^. Intuitively the 
phase-space directions where the wavepacket spreads are the unstable directions 
of the classical flow, whereas those along which they are squeezed are the stable 
ones ; those stable and unstable directions are encoded in the symplectic matrix 
Ft. 

All that follows will be true for a very general class of Hamiltonians (possibly 
time-dependent) : 

3m,M,K>0: {1 + z'^)-^'^/^\d2H{z,t)\ < K V|7| > m (3.2) 

uniformly for (t, z) G [— T, T] x Z 

such that the classical and quantum evolutions respectively (for the classical 
symbol and its Weyl quantization resp.) exist for t G [— T, T]. 
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It is well-known that the stability of the classical Hamiltonian evolution governed 
by H{z,t) is given by the following linear system : 

F = JMtF (3.3) 

where Mt is the 2n x 2n Hessian matrix of H at point Zt of the classical trajectory : 



is symmetric real, J is the symplectic matrix ()2.3jl and the initial datum is 

F(0) = 11 (3.5) 

Consider the purely quadratic Hamiltonian (time-dependent) : 

Ho{t) := ]^Z.MtZ (3.6) 

It induces a quantum evolution Uo(t,t') via 

th^^Uo{t,t') = Ho{t)Uo{t,t') (3.7) 

which is entirely explicit : 

Lemma 3.1 Let Ft be the 2n x 2n symplectic matrix solution of We 
note by R{Ft) the associated metaplectic operator, unitary in Tt. We have : 

Uo{t,0) = R{Ft) (3.8) 

and therefore by the chain rule : 

Uo{t,t') = R{Ft)R{Fi') (3.9) 

The proof of this result is classical and can be found in [Hllj . 

It encompasses the physical intuition that, for Hamiltonians purely quadratic, 
the quantum dynamics is exactly solvable in terms of the classical one. Namely the 
linear equ. ()3.3p is nothing but the classical Hamilton's equations for the quadratic 
Hamiltonian ()3.6|) . 

In fact R{Ft) decomposes itself into the product of two unitaries, (built from the 
symplectic matrix Ft) 
-one expressing the "squeezing" 
-one expressing the "rotation" 
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Lemma 3.2 We define the usual "creation and annihilation" operators of quantum 
mechanics as follows : 

a :— a' := a = (3.10) 

V2h V2h 



From Ft can be built two 2n x 2n matrices Et and Vt such that : 

R{Ft) = S{Et)R{t) (3.11) 

S{Et) := exp {^{a"^ .Eta"^ - a.E^a)^ (3.12) 

R{t) = exp (^^{alfta + a.F^a^)^ (3.13) 
Morally, in dimension n=l, R{t) has the simple form 

exp (f (Q^ + P^) 

which is simply a rotation by the real angle 7^. 

Let us consider the Taylor expansion up to order 2 of Hamiltonian H{z, t) around 
the point z — Zf {(lt,Pt) of the classical trajectory at time t : 

H^it) := H(zt,t) + (z- Zt).VH(zt,t) + ^{z - Zt).Mt{z - Zt) (3.14) 
By quantization it yields : 

H2{t) = H(zt,t)l+(Z - zt).VH(zt,t) + ^{Z - Zt).Mt{Z - Zt) (3.15) 

Let U2{t, s) be the quantum propagator for Hamiltonian H^it). We have : 
Proposition 3.3 

U^it, s) = e'^''-'^^'^f{zt)R{Ft)R{F;^)f{-Zs) (3.16) 

where 

5, := S,(z) - tPlpl (3.17) 

and St{z) = J ds{qs.ps — H{zs,s)) is the classical action along the trajectory z ^ Zf 
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The proof of this result can be found in jl2j . 

In fact this propagator which, being constucted via generators of the coher- 
ent/squeezed states acts in a simple manner on coherent states, and appears to be 
a good approximation of the full propagator UH{t,s) in the classical limit, when 
acting on coherent states : 

f/2(t,0) = e''^/''fizt)SiEt)R{t)fi-z) 
U2{tMz)=e'''"'f{zt)S{Et)R{tm 

where = \t''^^t 
The state 

^z,t) ■.= f{zt)SiEt)\0) (3.18) 
is simply a squeezed state centered in Zt, with a "squeezing" given by matrix E^. 

Theorem 3.4 Let H be an Hamiltonian satisfying the assumptions \3. jj)) and the 

existence of classical and quantum flows for t G [— T, T] . Then we have, uniformly 
for {t,z)e[-T,T]xZ : 

||t/H(t,0)(/?, -e'^*/^+^*$(2,t)|| <C^i{z,tY\t\Vhe{z,tf (3.19) 

P being a constant only depending on M and m, and 

fi{z,t) := Supo<s<t{l + \zs\) 

e{z,t) ■.= Supo<s<t{trF:Fsy/^ 
The proof of this result is contained in 

The estimate ()3.19|) contains the dependence in t, h, z of the semiclassical error 
term. One hopes that this error remains small when h ^ 0, provided that z belongs 
to some compact set of phase-space, and |t| is not too large. 

Typically 

e{z,t) ~ e*^ 

where A is some Lyapunov exponent that expresses the "classical instability" near 
the classical trajectory. The RHS of equ. ()3.19|) is therefore 0{h'^^'^) provided 

oA 

which is typically the Ehrenfest time, up to a factor 1/6 that is probably inessential. 
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Remark 3.5 (see fT^ ) Theorem (3.4) can be modified (and therefore also the state 
^{z,t)) to have an estimate in 



.7 = 1 ^ 



and therefore typically 0{h}^'^) with I integer as large as one wants. The squeezed state 
however now depends on I, and is typically a finite linear combination of wavepackets 
of the form : 

f{z,)R{F,)\^>,) 

where the are excited levels of the Harmonic Oscillator in dimension n. 

We shall now make use of the Mehlig- Wilkinson formula ()2.29|) . in order to 
calculate the dominant contribution to -R(t), the quantum overlap at time t . 

Of course this formula holds true as long as 1 is not an eigenvalue of F, but it 
can be generalized to this case also (see 

We now rewrite equ. ()3.19|) as : 

UH{tMz)=(^''''^T{zt)R{Ftm+e{tXz) (3.21) 
We shall consider the following quantity (return probability) : 

R{a,t):=\{^^,UH{t,Q)^a)\ (3.22) 

It can be rewritten as 

(f (a)(^o, f{at)R{Ft)^^) | + e{t, a, h) (3.23) 

The approximant ()3.23|) is used to obtain the following main result of this section 
which is a semiclassical approximation of R{a,t) ( the return probability) : Let us 
define the following 2n x 2n matrix 

K := {-J + tl){l-2iAy\j + il) (3.24) 

Theorem 3.6 There exists a constant ct, < q < 1 only depending upon F[t) 
such that R{a,t) has the following semiclassical approximation : 



Ct 



exp ( --^r(a - at).[l - (1 + iJ)(l - 2iA)-\l - iJ)]{ a — at 
An 



(3.25) 



The proof is rather technical, although not difficult, relying essentially on Gaussian 
integrations, and postponed to the Appendix. 

Remark 3.7 The fact that the exponential factor in \3. 2^) is a truly decreasing 
factors as a — at increases has been established in \T^ , using the properties of the 
matrix A ( essentially the fact that the eigenvalues of the matrix ^ — iA all lie on the 
circle \z — 1\ = 1, with z = excluded). 
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As regards the quantum revivals, we look at the expression ()3.25|1 . and where 
it equals 1 when time t goes on. Clearly we need at = a, and q = 1. The first 
condition is realized provided a lies on a periodic orbit 7 of the classical flow, and 
t = is the period of this orbit. The condition q = 1 is less transparent. It holds 
obviously if Ft = 1,, or more generally if Ft^^ = 11 cos ^ + J sin 6', 6 G [0, 2tt[, which 
is not automatically fulfilled along any periodic orbit, but can however occur for 
some t = pT^, p being a repetition number (integer). We thus have the following : 

Theorem 3.8 Let 'j be a periodic orbit of the classical flow corresponding to the 
(principal) symbol of H , wth period T^, and assume that z E '~f. Assume moreover 
that the stability matrix F^.^ takes the value 

Ft^ = tcose + J sine, ^G[0,27r[ (3.26) 

Then provided that < ^^\ogh~^ , we have that the overlap \Li{a,t) \ = \ {ipa,Utipa)\ 
taken at value obeys : 

R{a,t) = l-0{h') (3.27) 

This result expresses the semiclassical "almost recurrence" of the quantum state 
Uh(T^, 0)1-2) to 1^) for not too long periods of the orbit 7 on which the phase-space 
point z sits. 

Proof : 

It is not hard to check that if Ft^ has the form above, then the matrix N given by 
(6.26) equals 

and therefore | det = 1, whence q = | det Nt]^^^"^ = 1 (see (6.30)). Note that we 
have no general result ensuring that it is actually the case that ()3.26p holds true 
for any periodic orbit 7. The only result we know for sure is that it holds true with 
^ = if the classical flow is globally periodic (see [S]). However in the following 
section we also provide some explicit examples where the Stability Matrix can be 
shown explicitly to have the rotation form (|3.26|) . 

Remark 3.9 Actually we have established (see UJj) that |detA^| = 1, (and thus 
Ct = I) if and only if Ft^ is unitary, ( which is a generalization of the rotation form 
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4 Quantum Revivals for some time-dependent 
Hamiltonians 

As we have shown in Section 2, the time-dependent quadratic Hamiltonians play 
a special role in the understanding of Quantum Evolutions, since in particular they 
generate the Operators of the Metaplectic Group. They have the remarkable well- 
known property that the quantum dynamics is exactly solvable in terms of the 
classical dynamics. They are thus a good laboratory to investigate the behavior of 
the stabihty matrix Ft when it is propagated along a classical (closed) orbit. 

Thus in this section wc investigate the Return Probability of a coherent state 
for the quantum evolution Ut generated by a quantum Hamiltonian of the following 
form : 

Hit) = \z.S{t)Z (4.1) 

where Sit) is a real symmetric 2n x 2n matrix, and also particular non-quadratic 
Hamiltonians of the form 

^ ' 2 ■'^ ' 2 g2 
In the first case, the calculus of 

Kit) ■= (4.2) 

will be exact, and we make use of the result of Section 2 that the quantum propagator 
associated to Hamiltonian is exactly given by the metaplectic operator R{Ft), 
where Ft is the stability matrix solution of equ. (3.3) 

We clearly have : (by some abuse of notation we still denote by Ut{z) the Weyl 
symbol of the unitary quantum propagator Ut) 

iVa, UtVa) = j dzW^S^)Ut{z) = J dzW^,{z - a)Ut{z) (4.3) 

= ^-1 {W^>^{z)n{z)) {a) 
where J-' is the symplectic Fourier Transform, and 



Wliz) = {J'W^,){z) = {^o,T{z)<fo) ^ e-^'l'^ 



Therefore 



^{^WlX^) = 1 ..^^ X, ./o exp (-^z.(- - iA)z^ (4.4) 

^ |det(l-F)|-V2 ^\ 2h ^2 ' ) ^ ' 
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and thus 

F^{t) = \J^-\Wln){a)\ = |detiV,rV2exp|-^Ja.5iJa| (4.5) 

where ^ 

{I ~ ''^) = ~ ^^-^^ 

and 

Nt:^-i^{{l-iJ)Ft + l + iJ) (4.7) 

As shown in the Appendix (Lemma 6.2), | det Nt\ > 1, and Bt — whenever Ft — 1, 
so that in this case Fa{t) — 1. 

Analysis of a few more or less trivial examples : 

• H{t)^g{t) 



g being a T-periodic function of mean zero : 



T 



dtg{t) = 







Let G be the primitive of g that vanishes at t = 0. Then clearly the stability matrix 

F.= [\ J,,) (4.8) 
is T-periodic and obeys Ft — 1 together with Nt — —iJ so that 

Fair) = 1 

for any coherent state (pa-, namely tpa is recurrent along the periodic orbits of the 
classical flow (and this is even true without any phase). 

Note that the stability matrix F^ is always a rotation in this case : 

_ I cosG(t) sinG'(t) 
* " ' -sinG'(t) cosG{t) 



. H{t) = g{t)t-±^ 
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g being again T-periodic. Then, letting a — {q,p) G be any point in phase- 
space one gets by an exphcit calculus the return probability of a coherent state 



Utipa) = exp ( -^tG{t) - + q') sinG(t) - (/ + q') sin' (^^) ) (4.9) 

From this we deduce that if G{t) = dtg{t) = 2k7r, for A; e Z, then 



p2 ^ 



• H{t) = ^ + ^(Acos(wi) + ii)Q' 



The classical trajectories are the solutions of Mathieu equations, which are known 
to be either stable or unstable, depending on the zones in the parameter space 
(A,//) e E^. 

Assume (A, jj.) lie in some stability zone of the Mathieu equation. Then there 
exist a real parameter p together with an infinite sequence {cn\ne'z ^^at the 
stability matrix has the following form : 



h Cn cos[(2n + p)ujt/2\ Y.n sin[(2n + p)ujt/2] 



E„(2ri + p)cn sin[(2n + p)ut/2] ^ Cn(2n + p) cos[(2n + p)ujt/2] 

(4.10) 

where the constants C, D are determined by 

n 

D:^Y,^2n + p)cn 

n 

In general Ft is only queisiperiodic, but it is known that for the limiting curves 
A* — fk{^) of the stability zones, one has p — 2k, A; e N, so that in this case 

Fr = l, T = ^ 

UJ 

in which case we have exact revivals for any coherent wavepacket (pa- 
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Note that if the {cn}^^^, and p happen to satisfy 

(2 - up) ^ c„ = 2u; ^ 



nCr, 



then Ft reduces to a simple rotation matrix 



cos(7rp) sin(7rp) 
— cos(7rp) cos(7rp) 

Note that it has been estabhshed in Pj that beautiful recurrences can be obtained 
in this case in a much more general way if the initial wavepacket (coherent state) is 
chosen in an appropriate way, with a dispersion conveniently adjusted according to 
the classical solutions of Mathieu equations. Define 

LUD 

which is a positive real number, at least in the high frequency regime (and for 
parameters in the stability zones). Define 

V'(x) := (7rn)"i/^a-i/2 (^-a^ (4.12) 

be the reference normalized wavepacket. It is the ground state of a particular Ha- 
monic Oscillator of the form 

Its recurrences are defined through the overlap 

F{t) := {i;,Uti;) (4.14) 
It has been established in jH] that one has the exact result : 
Proposition 4.1 One has 

F{t) = e-'P^^'^Git) (4.15) 
where G is the T-periodic function defined as 

Therefore one has revivals at times T = ^ up to the secular phase exp (— 

|F(A;r)| = l, WkeZ 
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Remark 4.2 Note that the reference wavepacket ip is a kind of coherent (squeezed) 
state, localized around the fixed point x{t) = of the Mathieu equation, which is 
a periodic orbit reduced to a point. In the paper referenced above (JW), this point 
corresponds to the center of a quadrupole radiofrequency trap (Paul trap), which 
"traps" an ion in its center, classically as well as quantum mechanically. Ideally 
this holds for ever, but due to imperfections of the trap, and possible occurrence of 
several ions trapped together, there is an "escape time" which can be evaluated with 
respect to the various parameters of the ion and of the trap (see f77)/). 



H{t) = ^ + fit)^ + t- 



in dimension 1, g being a positive constant, and t \—>- f a T-periodic function. 
(This Hamiltonian has been considered in JU] as a model for interacting ions in 
a quadrupolar trap.) The important point is that the purely quadratic part of the 
Hamitonian governs the quantum dynamics. Let x{t) be a solution of the Hill's 
equation 

x{t) + f{t)x{t) = (4.17) 

with initial data 

x(0) = a, x(0) = - (4.18) 
a 

It can be shown easily that x{t) can be written as 

x{t) := e"+'^ (4.19) 
where t ^ u,6 are real functions, and u is T-periodic, and that 

e = e~2« 

provided we are in a stability region of equ. ()4.17|) (which are generalization of the 
stability zones of the Mathieu equation considered above). Then it can be estab- 
hshed : (see [lU]) 



Proposition 4.3 Define 
and for any n G N 



:= ^VT+8? (4.20) 



E„ = 2n + a + 1 (4.21) 



^n(x) = ( ^, l""-^.. x^^^'^e-^^'^Llix') (4.22) 
V[a + n + l) J 
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L° being the Laguerre polynomials. Then a solution of the time- dependent Schrddinger 
equation for H{t) is : 



^.(x, t) = exp (^-tOE^ + ^i±!!^^ ^„(xe-") (4.23) 

Thus, since u is T-periodic, and 'u(O) = due to our choice of initial data ()4.18j) . this 
result expresses the recurrence, up to a phase, of the corresponding set of wavepack- 
ets. 

Proof of Proposition : 
It is easy to check by direct computation that the unitary evolution operator corre- 
sponding to Hf is of the form t/tf/g, with 

JJ^ ^ ^iuQy2^iu{P.Q+Q-P)/2^-ieHo (4.24) 

where we make use of the following differential equation satisfied by u : 

M + M2-e-^" + / = (4.25) 



n 

5 Quantum Fidelity in the linear-response semi- 
classical regime 

Consider the following quantum Hamiltonians Hq and H := Hq + W, where V 
is some perturbation, and A is a (small) coupling constant. We want to compare 
the quantum evolutions Uuoif) under Hq, and Unit) under H respectively. Starting 
with a reference normalized state ip, the quantum fidelity amplitude (or Loschmidt 
Echo), is defined as : 

Fit):=\{UHom,UHm)\' 

Of course if t = this quantity is equal to 1, and the fidelity is assumed to decrease 
as time increases. The point is to obtain suitable regimes where something can be 
said about this decrease. 

If we choose the reference state ip to be an eigenstate of either Hq or H, then 
clearly the fidelity reduces to the return probability for either H, or Hq. In 
this sense the studies of Return Probability and Quantum Fidelity are linked. 



20 



As an entertainment exercise, we study the quantity F{t) when i/j is taken as an 
eigenstate ipj of the unperturbed Hamihonian Hq corresponding to discrete eigen- 
value Ej, and recall an old (exact) estimate on the short time fidelity amplitude 
squared in that case. 

F^it) = \{i;„UHm,)\' 
which means that we only have to consider the "return probability" in this case. 

The result we recall here is known as the "Mandelstam-Tamm inequality" |2H]- 
(see also j2ZI) : Let 

E:={i),H^lj) (5.1) 



AE^ := {ij,{H ~ Efij) (5.2) 

Proposition 5.1 Assume E, AE > are defined by equ. h5. 1\) and h5. g)) . respec- 
tively. Then the return probability of a state ip with respect to the quantum evolution 
Unit) generated by the Hamiltonian H, \ {ip, UH{t)ip)\^ remains, as long as time stays 
bounded above by 'Kh/2AE, bounded from below by cos^(^^) .■ 

|(^,MtWr>cos2(^) (5.3) 



One of the most studied regime is the so-called "linear response" regime, which 
is a perturbative one, where the Dyson series of Unit) is truncated up to second 
order expansion terms. Let us denote by V{t) the Heisenberg observable : 

V{t):=UH,{-t)VUH,{t) (5.4) 

and by V{t) the "mean" of this operator on the interval [0, t] : 

V{t) := [ V{s)ds (5.5) 

One defines : 

/^^(t) ;= 1 + V{t)^) - ^ £ ds{^, V{s)Vis)^) (5.6) 

Since we are mainly interested by the absolute value (squared) of this quantity (to 
obtain the "quantum fidelity") , we now consider F^^(t) to be the modulus squared 
of f^^it) truncated to second order in A : 



For such a truncated Linear-Response Fidelity, we shall perform the small Planck 
constant limit, namely the semiclassical regime under suitable assumptions on Hq : 
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We shall now take as particular reference states the eigenstates of Hq whose 
eigenvalues lie in some neighborhood of a given classical energy E. Let us consider 
h- dependent energy intervals in which the spectrum of Hq is pure point : 

lih) := [aih),m] 

where a{h) < E < [3{h) with \mi}i_^Q{[3{h) — a{K)) = 0, P{h) — a{h) > Ch for some 
C > 0, 

and denote : 

m := {j : G m} 

Let ipj be an eigenstate of Hq for some j G A(/i). For such reference state, the 
expression in (5.7) equals : 

F^^'it) = 1 - ^ {2^ j\s j^^ ds'{^,,VUH,{-s')VUH,{s)ij,) - V(t)V^,f ^ 

(5.8) 

We make use of the following result : 

Proposition 5.2 Let A he the Weyl quantization of a classical symbol A, and let 
us consider the Wigner transform associated with a given state ip ^Ti : 

W^{q,p) := h-^ [ dyi,{q + |)V^(g - he'^-vl^ = /i"" [ dz' {^P,fiz')^P)e''^^'^''^^^ 

(5.9) 

Then we have : 

(V',iV)= / dzW^{z)A{z) (5.10) 
Jz 

Let V{z) (resp. Ag) be the principal symbol of V (resp. VV{s)). Consider the clas- 
sical flow induced by the classical symbol of Hq. Then we have the following 
semiclassical result (Egorov theorem) : 

Proposition 5.3 As h ^ {] we have : 

AXz)^V{z)Vo^^^^X^) (5.11) 

uniformly for (s,z) in any compact set o/M x Z . 
This result goes back to Egorov [TH] (see also P) 

Now the Wigner function of (almost any ) ipj, for j G A{h) can also be shown 
to have very interesting semiclassical limit properties, for ergodic classical flows, 
namely to converge towards the microcanonical measure on the energy surface Hq = 
E as h in a weak sense (known as Schnirelman Theorem, see ^Hl) : 

Ho{q-,p) being the principal symbol of Hq., we denote by the energy surface 

:= {(g,p) e , HQ{q,p) = E] (5.12) 
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and by daE the microcanonical measure normalized to unity : 

where dT,E is the Euchdean measure on Eg. 
We make the following assumption : 

(H) The dynamical system {T^e, daE, is ergodic, which means : 
for any continuous function a on E^;, we have for almost every 2; G : 

lim T"^ / ao^^jj{z)dt= ! a{z)daE{z) (5.14) 
Jo Jt.e 



Proposition 5.4 Under rather technical assumptions (see and hypothesis (H), 
for any h > 0, there exists M{h) C K{h) depending only on Hamiltonian Hq such 
that : 

and 



lim Ul).,Ai).) = / Aiz)daE(z) (5.16) 



E 



The proof of this proposition has been given in |T^ . 

We shall employ all these results to give the semiclassical behavior of the linear- 
response quantum fidelity. 

We define the mean and the autocorrelation function of V (on the energy shell) 
as follows : 

Ve:= I daE{z)V{z) (5.17) 



CvAt)-= f daE{z){V{z)Vo¥^^{z)-Vi) (5.18) 

Theorem 5.5 Under assumption (H) we have the following result : 

hm ^''Z\T^'^ =1-4 fds'CvAs') (5.19) 

^-0 ^k{h) Jo Jo 

The proof of this result follows easily from Propositions 5.3-5.4 applied to equ. 
(5.7). 
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Remark 5.6 This result has first appeared in with a more heuristic approach. 
Note that useful estimate of the second term in the RHS of i5.19\} could be obtained 
under better knowledge on the possible "mixing properties" of the system (and there- 
fore on the decay of correlations Cv,E{t) as t becomes large). 

Usually physicists consider only the two first terms, claiming that it can be "ex- 
ponentiated" to yield an "exponential-type" decay of fidelity (in the linear-response 
semiclassical regime) of the form exp{—CX'^t/h?), provided the "mixing" is strong 
enough. C is a constant given by the mixing property of the classical flow. The ar- 
gument goes as follows : 

If the decay of correlation function is integrable, (for example if the mixing is ex- 
ponential) then the integral over s' in equ. (5.19) is bounded uniformly in s, and 
therefore the RHS of (5.19) is bounded above by 

They now "exponentiate" this estimate to yield the "large-time" behavior of Fj"^{t). 
We are not able to go further now in that direction by our rigorous approach. 

6 Appendix 

We give now the proof of Theorem 3.6 : Using repeatedly the following property 
of the Weyl-Heisenberg operator T{z) : 

Lemma 6.1 i) 

f{z)f{z') := e-^-(^'^')/2^f + z') 

ii) 

{^q,T{z)<^q) = exp(- — ) 
we get easily that the RHS of equ. ()3.23|) equals (up to the error term e(t, a, h)) : 

|(<^o, f{at-a)R{Ft)^o)\ (6.20) 

Now using the Mehlig- Wilkinson formula (2.25) for R{Ft), we get, provided 1 is 
not eigenvalue of Ft : 



\{ifo,fiat-a)R{Ft)ipo)\ = | det(l - F^r^/^ 



detd-K)!-!/^ 



(6.21) 

du{ipo,f{u + at- a)<^o)e^-^"/2/ig-i-K-a,«)/2a 

We use the simplifying notation : 

Zt := at — a (6.22) 
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Now, forgetting the factors h, that we shall reintroduce later by homogeneity, we 
let that equ. ()fi.2()j) equals 



1 



det(l-Ft)|V2 
1 



det(ll-Ft)|i/2 



duexp { -u.Jzt + 1—^- ) {fo,T{u + Zt)(po) 



. iu.Jzt 1 , ,0 ^ 
du exp ( — -{u + Zt) + -u.Au 



(6.23) 



= |det(l-Fj|-i/2g-i 
It is easy to check that the 2n x 2n matrix ^ — iA obeys 



^{--lA) > 

so that by the formula of Fourier Transform of Gaussian functions, this yields : 



det(- - iA) 
^2 ^ 



-1/2 



det(ll-F,)|-V2 



exp (^-^ - w.{l - 2iA)-'^iu^ 



(6.24) 



where w denotes the following complex vector in C 



2n 



w := ^{J + il)zt 



Thus ()6.24|) can be rewrtitten as 



det{--tA){l~Ft] 



-1/2 






exp 



1 



where we have defined K as in equ. (3.24) 

Recollecting the various prefactors, and reintroducing the parameter h, we finally 
get as a semiclassical approximant of R{a, t) the following expression : 



R{a,t) 





-1/2 




iA) 




exp ^ 



1 



(6.25) 



We have the following Lemma : 
Lemma 6.2 i) We have : 

1 



25 



where N is a 2n x 2n matrix such that \ det > 1 

a) The matrix K is such that : 

I exp(2;.i^z)| < exp(z^) 

for all 2 e M^". 
Proof : 

I - = i-iJiF -1) + F + 1){F- 1)-^ = N{F - 1)-^ where 

N ■= -^-^{{l-iJ)F +l + iJ) (6.26) 

which has the 4-block form : 

N = -—( ^ + ^ ^' + '^ 
2 V ^(^ - 1) *(^' - 

where Z = a — ic, Z' = b — id and a, 6, c, d are the n x n block matrices of 
which F is composed : 




But using the symplecticity of F, the resuh that | det > 1 can be easily checked. 

ii) is established in jHj. It uses the fact that w.w = and that ^ — iA has all 
eigenvalues located on the circle \z — 1| = 1, with origin excluded. □ 

Then we conclude that the prefactor in ()6.25|) is q := |detA^|~^/^ < 1, which 
completes the proof of Theorem 3.6. □ 



7 CONCLUDING REMARKS 

Our work is a first mathematical approach of estimates for the Quantum Return 
Probabilities (and Quantum Revivals), and of the Quantum Fidelity for a rather 
large class of Hamiltonians. It also contains a rigorous proof of the Mehlig- Wilkinson 
formula for the Metaplectic Transformations, which appear to be an useful tool, 
when combined with the coherent states framework, to understand the semiclassical 
behavior of the Return Probability. Then we are able to describe under which condi- 
tions on the linearized flow we can have Quantum Revivals for wavefunctions being 
Coherent States located on a periodic orbit of the classical flow. We stress that 
a similar tratment can be performed for a semiclassical analysis of the Quantum 



(6.27) 
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Fidelity. This will be developed elsewhere [Ti] . 

Here we have concentrated our attention on the Linear Response Regime for the 
Quantum Fidelity. It has been already studied in the physics literature by several 
authors ([H1,|S2I, [SS], IS])- Furthermore the link of the Fidelity Decay with the 
question of Decoherence has been put forward ([IS], [T7j, [201, |2S]; |2E]; [SI])- 
It is an interesting topic that deserves further investigation, but that we have not 
worked upon. 

Since no complete review is up to now available on Quantum Fidelity, contrarily 
to the Quantum Revivals problem j^, we have provided in the bibliography, for 
the readers' convenience, a rather long (however not exhaustive!) list of references 
on the physics literature on Quantum Fidelity. It is not the purpose here to discuss 
them in detail. To my opinion, one of the prominent points in these approaches 
is that they interestingly distinguish several regimes (Fermi-Golden-Rule, Semiclas- 
sical. Linear Response) where significantly different time behavior of the fidelity 
manifest themselves ; moreover they study the link between these behaviors and the 
regular/chaotic nature of the underlying classical motion. We think that these topics 
open new research lines in mathematical physics. 

Let us stress finally that the results developed in this paper, regarding Quantum 
Revivals as well as Quantum Fidelity, show a delicate interplay between time t, 
Planck constant, and the size of the perturbation A. In particular in the 
Linear Response Regime, the parameter that governs the fidelity decay is the com- 
bination tX/h. Regarding the Semiclassical Regime, we do not seem up to now to 
be able to go beyond the so-called "Ehrenfest time" C\ log^|. 
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